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B decays into tt-kK and KKK, where the -nn and KK pairs interact in isospin zero 5-wave, are studied in the 
7T7T effective mass range from threshold to 1.2 GeV. The interplay of strong and weak decay amplitudes is analyzed 
using an unitary mr and KK coupled channel model. Final state interactions are described in terms of four scalar 
form factors constrained by unitarity and chiral perturbation theory. Branching ratios for the B — » /o(980)A 
decay, calculated in the factorization approximation with some QCD corrections, are too low as compared to recent 
data. In order to improve agreement with experiment, we introduce long-distance contributions called charming 
penguins. Effective mass distributions, branching ratios and asymmetries are compared with the existing data 
from BaBar and Belle collaborations. A particularly large negative asymmetry in charged B decays is predicted 
for one set of the charming penguin amplitudes. 



1. Introduction 

Recent experimental results from B factories 
indicate that charmless hadronic three-body de- 
cays are more frequent than two-body ones [1]. 
Moreover one observes on Dalitz plots a defini- 
tive surplus of events at relatively small effective 
masses. This is a signal of especially strong inter- 
actions between hadrons at not too high relative 
energies. Many resonances are explicitly visible 
but in general the interference pattern is quite 
complicated. Knowledge of these final state in- 
teractions is important to obtain a precise de- 
termination of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. Weak decay observablcs 
give information on hadron-hadron interactions 
and internal quark or hadronic structure of the 
produced particles. 

Prominent maxima in the tt + it~ spectra are 
observed in the B — > ir + ir~K decays in the 
/o(980) region [2-9]. The B+ -> f (980)K+ 
and B° — > /o (980) AT branching ratios arc rel- 
atively large, and of the order of 10~ 5 . Direct 
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and time-dependent CP- violating asymmetries 
are also measured. The first Belle result on the 
B + — > /o(980)A"+ branching ratio [2] has moti- 
vated the study of Chen [10]. In the perturbative 
QCD approach Chen finds that the non-strange 
content of the /o(980) can be important. Accord- 
ing to Cheng and Yang [11], sublcading correc- 
tions due to intrinsic gluon effects inside B meson 
may enhance the decay rate of B — > / (980)A'. 
B decays into scalar-pseudoscalar or scalar- vector 
particles have been studied by Minkowski and 
Ochs with a special emphasis on the presence of 
the lightest glucball [12]. The 7r + 7r~ mass spec- 
trum in B — > A'7T7r decays, reported by Belle 
in 2003, is reproduced by a model amplitude of 
the coherent sum of /o(980), /o(1500) and a very 
broad glueball as a background. 

In the present letter we study the B decays 
into ttttK and KKK. We restrict ourselves to 
the case where the produced tttt or KK pairs in- 
teract in isospin zero S'-wave from the 7T7t thresh- 
old to about 1.2 GeV. One expects the tttt isospin 
two 5-wave contribution to be small since the up- 
per limit of the branching fraction for the B + — > 
tt+tt+A"- decay is less than 1.8 x 10~ 6 [13]. Us- 
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ing the KK/nrj branching ratio of ao(980) [1] and 
the upper limit of 2.5 x 1CP 6 [14] for the branch- 
ing ratio of B+ -> ag(980)AT+, ag(980) -► tt ??, 
one can estimate the branching fraction B{B + — > 
a%(980)K+, ag(980) -> K+R-) to be smaller 
than 1 x 1CT 6 . This indicates that the KK 
isospin one 5*- wave amplitude is suppressed in the 
B ± -> K+K-K ± decays. 

Two-pion iS-wave rescattering effects have been 
recently considered by Gardner and Meifiner [15]. 
They study the effect of the /o(600) (or a) reso- 
nance on the B° decay into 7r + 7r~7r° in the range 
where the p(770)ir channel dominates. The <77r 
channel can play a role in the determination of the 
CKM angle a from the B° — ► pir decays. Gard- 
ner and Meifiner describe the broad /o(600) in- 
troducing a scalar form factor constrained by the 
chiral dynamics of low-energy meson-meson in- 
teractions [16]. This scalar form factor is used in- 
stead of the commonly applied Breit-Wigner form 
to improve the description of the broad a and the 
understanding of the B — > pir decays. 

We extend the approach of ref. [15] to the 
/o(980) resonance. The four strange and non- 
strange 7T7T and KK scalar form factors are con- 
strained by chiral perturbation theory as devel- 
oped by Meifiner and Oiler [16]. Our final state 
interaction is, however, different from that of [15] 
and [16]. Here we consider the unitary ttt: and 
KK coupled channel model of [17]. 

First the B — ► (tttt)s— wave 

K and B -> 

(KK)s-wave K decay amplitudes are calculated 
within the naive factorization approximation [18, 
19]. Penguin amplitudes interfere destructively 
which leads to much too small B — > /o(980)A' 
branching ratios. Then we consider some 
QCD factorization corrections [20] calculated by 
de Groot, Cottingham and Whittingham [21]. 
These corrections are not sufficient to obtain 
agreement with experiment. Further contribu- 
tions are needed. Here we include the long- 
distance contributions which have been consid- 
ered in [21] to improve their fit to hadronic charm- 
less strange and non-strange two-body S-decay 
data. These amplitudes, called charming pen- 
guin terms, originate from enhanced charm quark 
loops [22]. They could, for instance, correspond 



to weak decays of B to intermediate DgD^ 
states followed by transitions to /o(980)i^ final 
states via cc annihilations. Their addition allows 
us to obtain a good agreement with the measured 
B — > /o(980)if branching fractions. 

In sect. 2 we describe our weak decay ampli- 
tudes supplemented by the scalar form factors. 
Our model for the final state interactions is given 
in sect. 3. Results of calculations and compari- 
son with available data are presented in sect. 4. 
In sect. 5 we give some conclusions and final re- 
marks. 



2. Amplitudes for the B 
KKK decays 



irirK and B 



We shall write the model amplitudes for the 
following decays: B ± — > (xTr)sK^, B ± — > 
{KK) S K±, B° (irir) s K°, S ->(KK) s K°, 
B° -> (TTTT^i^and B° -> (KK) S K°. Here by 
(irir)s and (KK)s we mean tt + it~ or 7r°7r° and 
K + K~ or K°K pairs in isospin zero S'-wave. 

The possible quark line diagrams for the 
B~ decay, together with the final state mesons, 
are shown in fig. 1. For the B° decay there are 
only two types of penguin diagrams similar to 
those shown in figs. 1(b) and 1(c). The tree dia- 
gram of fig. 1(a) is absent. The uu or ss transi- 
tions into 7T7T or KK states, shown in fig. 1, are 
described by four scalar form factors. 

In the approximation used in our approach the 
effective weak Hamiltonian H is replaced by the 
sum of products of factorized currents [19]. We 
introduce some of the QCD factorization correc- 
tions and the charming penguin amplitudes con- 
sidered in [21]. Then the B~ -> {tt+tt~) s K' 
decay amplitude is 

((Tr+7r-)sK-\H\B-) 

= %yf {x [P{m^)U + C{m^)] rnm™)+ 

+ [Q(m ww )V + X C{m K )}Tl*(m^)}, (1) 

where Gf is the Fermi coupling constant and x 
is a constant which will be estimated from the 
properties of the /o(980) decay. The functions 
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Figure 1. Quark line diagrams for the B~~ de- 
cay: (a) tree diagram, (b) and (c) penguin dia- 
grams. The spring-like lines represent gluon ex- 
change and the dashed ones the irir or KK isospin 
zero ^-wave pairs. 



and F^~" K (m^ v ), respectively. The masses 
of B meson, kaon, pion, b-quark, strange-, 
down- and up-quarks arc denoted by Mb, Mr, 
m w , m b , m s , rrid and m u , respectively. In 
eq. (2) Jk is the kaon decay constant. In cq. (3) 
Bq is related to the vacuum quark condensate: 
Bo = -{0\qq\0)/f%, U being the pion decay con- 
stant equal to 92.4 MeV. We use the formula 
Bo = mj/(2ra), where m is the average mass of 
the light quarks u and d. We put m = 5 MeV and 
following ref. [19] we take m s = 0.122 GeV and 
m b = 4.88 GeV. The functions U and V in (1) 
depend on the combinations of the coefficients at 
[18-21] and on the products of the CKM matrix 
elements X u = V ub V* s and A t = V tb V t * s : 

U= X u [oi + a\- a% + {a% - dg) r] + A t (a c 6 r~ a%) , (4) 



V = X u (a c 6 - <) + X t a c 6 , 



(5) 



where the chiral factor r = 2M^/[{m b +m u ){m s + 
m u )]. In the numerical calculations we set m u = 
fh and the coefficients a\, d%, a%, ag and a%, eval- 
uated at the scale \i = 2.1 GeV, are taken from 
table III of [21]. These coefficients take into ac- 
count some QCD factorization corrections. We 
do not include small corrections coming from hard 
gluon exchanges with spectator-quark, from anni- 
hilation terms and from electroweak penguin di- 
agrams. The charming penguin contribution can 
be parametrized as 



C(m) 



(M B - 



)f*F v (A u P 1 G/M +A t Pi) , (6) 



r™(m nn ) and rf(m WT ) are the non-strange and 
strange pion scalar form factors depending on the 
effective pion-pion mass m^. Furthermore the 
functions P{m^ 7I ) and Q(m vn ) defined as 

P{m^) = f K {M% - m^)F^^{Ml), (2) 
2^50 



m b -m s 
are proportional to the B 



{Ml-Ml)F B ^ K {ml«), (3) 



B —> K transition form factors, F Q 



(tttt)s and 
(Mi) 



B — >(Tnr) s 



where m is m T7r or Mk, F„ is the B — > n tran- 
sition form factor calculated at the zero m^. limit 
and F™ 7 , Pi are complex parameters. Determi- 
nation of these parameters has been done in [21— 
24] by fitting some charmless two-body i?-decay 
data. 

The neutral B — > (iT + ii~)sK'g amplitude is 
similar to the charged one: 



V2 



with a± = and m u — > m^. 
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The amplitudes for _B-dccays into three kaons 
read: 



{(K+K-) S K-\H\B-) 

x [rr(rn K K)+T?(m K x)] 
+ [Q{m K K)V + X C{m K )} 

x m*(m K K)+r s 2 *(m Kl ?)}} (8) 

and 

((g ^- )3g o| g |B°> = <g^M1gE (9) 



with a\ — 0, F 2 ' s (m K -w) 



V2 
and m. 



m d . 



Here m K -^ is the K + K 



effective mass and : 



is the effective mass with the second possible 
K + K~ combination. 



T 2( m KK) and T 2( m KK) 

are the non-strange and strange kaon scalar form 
factors. Decay amplitudes with (tt°tt°)s and 
(K°K )s final states are the same as those with 
(ir + ir~)s and (K + K~)s pairs, respectively. 

If one uses the ai coefficients [19] a| = a\ = 04 
and ag = = a6, then our formula (1) with 
C(m) = has the same algebraic structure as 
the B~ — > <77r~ amplitude of rcf. [15] (sec their 
eq. (25)). The above equalities for are valid in 
naive factorization but not in QCD factorization. 
A particular feature of the two penguin contri- 
butions to the b — > s transition is the near can- 
cellation of these two terms in eq. (4) due to 
a% f» cig and r ~ 1. This has been pointed out 
by Chernyak in his estimation of the scalar pro- 
duction in B decays [25] and also by Gardner and 
MeiBner [15]. 

Replacing the A„ and At values by their com- 
plex conjugate values A* and A^ in eqs. (1), 
(7), (8) and (9) gives the amplitudes for the 



B+ -> (ir+Tr-) s K+, B (] 
{K+K-) S K+ and B° - 



■)sK° s , B+ 



{K+K-) S K% decays. 



3. Final state interactions 

In the B decays considered above, one should 
include in the final states the tttt — ► tttt or the 
KK — > KK rescattering and the tttt — > KK 



or the KK — > tttt transitions. The (tttt)s and 
(KK)s pairs are formed from the uu, dd and ss 
pairs. The four scalar form factors appearing in 
eqs. (1) and (8) are defined [16] as 



r?(m)> 

Tf(m)> 
T s 2 (m)j 



1 



V2S \(0\nn\KK) 



^2B \(0\ss\KK) 



(10) 



where to is the effective tttt or KK mass, rin = 
(uu + dd) I y/2 and |0 > denotes the vacuum state. 
The final state interactions, which satisfy the uni- 
tarity constraints, are incorporated in the follow- 
ing formulae: 



(11) 



where and \kj) represent the wave functions 
of two mesons in the momentum space and the 
indices i, j = 1,2 refer to the 7T7T and KK chan- 
nels, respectively. The center of mass channel 
momenta are given by k\ = y m2 /4 — m 2 and 
fc 2 = y to 2 /4 — m 2 K . The matrix T is the two- 
body scattering matrix. Here we use the solution 
A of the tttt and A' A coupled channel model [17]. 
As we restrict ourselves to m T7r < 1.2 GeV, 
the third effective (27r)(27r) coupled channel con- 
sidered in this model, with a threshold around 
1.4 GeV, has a small effect. The functions Gi(m) 
are the free Green's functions defined in [17] and 
A™' s (to) are the production functions responsi- 
ble for the initial formation of the meson pairs 
prior to rescattering. The production functions 
have been derived by MeiBner and Oiler in the 
one- loop approximation of the chiral perturbation 
theory [16]. Using their eqs. (37), (38), (41) and 
(44) one obtains 



ATM = 
i?"(m) = 
R\{m) = 
R s 2 (m) = 



0.566 + 0.414m 2 , 
-0.322 + 0.527 to 2 , 
-0.036 + 0.353 to 2 , 

0.071 + 0.338 to 2 , 



(12) 



where to is in GeV. 
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If one considers only the on-shcll contribu- 
tions in cq. (11) then the scalar form factors can 
be written in terms of the phase shifts 5^(m), 
$Ki((jn) and of the inelasticities 77(771): 



ry (m) 
i 

~~ 2 



R^ s {m) (l + ri{m)e 2iS ^ m A - iR%' s (m) 



x Jp-Jl- ^( TO ) e *[«~("0+W"0] 
V ki 



, (13) 



ry (to) 



R^ s (m) (l+r;(TO)e 2l ' 5 ^ (m) 



iR™ ,s (m) 
(14) 



— J\ - rj 2 (m)e i[57 ' 7 ' W+W m )] 
k 2 



Below the KK threshold 77(777) = 1 and 

ly s *(m) = i?"' s (TO)cos^(m) e lS ™ {m) . (15) 

Particularly interesting is the m^r range where 
the phase shifts are close to 180°. Then one 
expects a maximum in |r™' s |. As we shall see in 
the next section, this is the case for the produc- 
tion of the /o(980) resonance. Note that the T™' s 
is zero when ~ ir/2. 

The off-shell contributions to the form factors 
are very much model dependent and are not con- 
sidered in the following calculations. 

4. Results 

The amplitudes for the B — ► (inr)sK decays 
considered in sect. 2 depend only on the effective 
mass m nn . Integrating on the Dalitz plot over the 
kinematically allowed range of m^K-, one obtains 
the differential B — > ttttK decay distribution 



dT 



h pk 



dm** 4A/|(2tt) 



x \M(B^> (mr)sK)\ 2 , (16) 



where ki and px = y/Ej^m^) — m 2 K are the 
pion and kaon momenta in the rnr center of mass 
system, E K (m V7T ) = \{M% - m 2 ^ - m 2 K )/m V7T 
being the corresponding kaon energy. In (16) M. 



denotes the decay amplitude given by eq. (1) or 
(7). Dividing dT/dm^ by the appropriate B + 
or B° total width T b one obtains the differential 
branching ratio dB/dm nn . 

Before presenting our results we fix the con- 
stants which appear in the formulae for the decay 
amplitudes (1), (7), (8) and (9). The masses of 
pions, kaons, i?-mesons and their life times, the 
values of the Fermi coupling constant Gf and the 
kaon decay constant fx = 0.1598 GeV are taken 
from [1]. For the kaon mass we use the average 
of the charged and neutral ones. Following the 
results obtained in [26] and applied in [15] we use 
F B ^ (7r7r)s (M^)=0.46 for the B {nir) s transi- 
tion form factor. 

The B — > K transition form factor is ap- 
proximated by a constant equal to 0.39 which is 
close to the number F^~* K (0) = 0.379 quoted 
by Bauer, Stech and Wirbel in table 14 of [18]. 
This approximation is justified since we consider 
a relatively narrow range of the effective masses 
m n7r . These masses are much smaller than the 
mass of the heavy excited B meson used in po- 
lar models of the transition form factor. So we 
fix Fq ^ k (m 2 ^) « Fq ^ k (m 2 o ) where to/ is the 
/o(980) mass. The value 0.379 quoted in [18] as 
well as a more recent value of 0.33, obtained by 
Ball and Zwicky in [27], are within the limits 
given by Beneke and Neubert in table 1 of [28] 
for F^ K (0) = 0.34 ± 0.05. 

Using the Wolfenstein representation [29], the 
CKM matrix elements are written in a form accu- 
rate to the level of A 6 [30] . The values of the pa- 
rameters, taken from the CKMfittcr Group, are: 
A = Vus = 0.2265, p = 0.189, rj = 0.358 and 
A = 0.801 [31]. 

The constant \ t> e fitted to the exper- 
imental branching ratio of the decay B + — > 
/o(980)A' + . Once fixed it will be used to make 
absolute model predictions in the whole TO W7r 
range studied here for the i? ± — > n + n^K ± re- 
action and for other decays like B° — ► ir + Tr~Kg, 



B+ K+K-K+ and B° K + I<-K° s . The 
value of x is, however, not arbitrary. It can be es- 
timated from the following considerations. First 
we shall concentrate ourselves on the m n7r range 
close to the relatively narrow resonance /o(980) 
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clearly visible in the B decays into tt + tt~ K. The 
/o(980) decays mainly into tttt. The coupling con- 
stant of Jo (980) to the tttt pair can be approxi- 
mated by 



= m/ 



'8 7rr(/o 



fci(m /o 



(17) 



(see, for example, eq. (36) of [32]). Here r(/o — ► 
7T7r) is the 7T7T partial width of the /o(980) and 
ki (m f ) is the pion momentum in the /o (980) rest 
frame. The two scalar form factors ^(mmr) and 
r* (m W7r ) are strongly peaked at the / (980) mass; 
one can show, however, using eqs. (12) and (13) 
that |r"(m/ )| 3> |r?(m/ )|. Therefore one can 
write the approximate relation 



X 



9 fan-" 



m fo r tot (fo) |ry(m/ )| 



(18) 



(see, for instance, eq. (35) of [15] for the /o(600) 
case). Here r tot (/ ) is the total /o(980) width. In 
ref. [32] r tot (/ )=(71± 14) MeV. For m fo = 0.98 
GeV we obtain from eq. (13) |Pf (m /o )| w 0.96. If 
Ttot(fo) « r(/ -> tttt) = 60 MeV then eqs. (17) 
and (18) give X ~ 30 GeV -1 . 

In this paper we do not attempt to make our 
own fits to data by adjusting the constants P^ IM 
and Pi in the function C(m) of (6). We present 
the results obtained with the charming penguin 
amplitudes determined in [21] and [24]. The the- 
oretical curves shown in figs. 2 to 5 correspond to 
the first set of amplitudes. 

4.1. — y 7T + 7v~K :iz decays 

In figs. 2 and 3 we show a comparison of the 
tttt effective mass distributions for the B ± — ► 
TT + TT~K ziz decays calculated in our model with the 
results obtained by BaBar [5] and Belle [7]. 

The branching fraction B(B ± — > 
/o(980)K±, / (980) TT+TT-) = (9.2±1.2±|i)x 
10~ 6 obtained by BaBar [5] can be reproduced in 
our model for the value of x — 35.0 GeV -1 . 
This value is close to our estimation given 
above. The Belle Collaboration has reported 
a slightly smaller value of the branching ratio 
(7.55 ± 1.24 ± 0.69to'.9e) x 10 ~ 6 I 7 !- Both valucs 




0.95 1 1.05 



Figure 2. 



Effective tt 



+ n 



mass distribution in 



7r + 7r K ± decays. The BaBar data are 



taken from [5]. 
model. 



The solid line results from our 



are compatible within their error bars. The av- 
erage value given by the Heavy Flavor Averaging 
Group (HFAG) is equal to (8.49^;^) x 10~ 6 [33]. 
For this value the constant \ is 33.5 GeV -1 when 
we use the charming penguin amplitudes of [21]. 

In fig. 2 the normalization of the theoretical 
curve to the data in the /o(980) range is based 
on the total number of events seen by BaBar for 
from 0.9 to 1.1 GeV, multiplied by a correc- 
tion factor of 0.92 being a ratio of the branching 
ratios 8.49 x 10~ 6 (HFAG's value) and 9.2 x 10~ 6 
(BaBar's value). 

In fig. 3 the solid curve is normalized at 
jtItttt = 976 MeV. This corresponds to the max- 
imum of the background subtracted mass distri- 
bution. We calculate it from fig. 9e of [7] as 
133 events/50 MeV. This number is obtained by 
taking into account two factors. The first fac- 
tor, equal to 92%, follows from the fraction of the 
/o(980) components in the full spectrum result- 
ing from the solution 1 of the Belle model, called 
Ktttt — C [7]. The second factor, equal to 1.125, 
comes from the ratio of the above mentioned 
branching ratios: 8.49 x 10~ 6 and 7.55 x 10~ 6 . 
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Figure 3. Comparison of the Belle [7] ir~ 
tive mass distribution for the B ± — > 7r + 7 
cays with our model. The dashed line corresponds 
to the B + — ► 7r + 7r _ A' + decays, the dotted-dashed 
line is for the £>~ — > 7r + 7r~A'~ decays and the 
solid line is the average for the B + and B~ de- 
cays. 



One can see from fig. 2 that our model de- 
scribes quite well the tt + tt~ spectrum measured 
by BaBar in vicinity of /o(980). The model de- 
picts also a very pronounced maximum seen by 
Belle near 1 GeV (fig. 3). This maximum is at- 
tributed to the / (980) resonance. At lower nn 
masses near 500 MeV one can notice a broad the- 
oretical maximum which we can relate to the a or 
/o(600) meson. Experimental data in this mass 
range are not in disagreement with this feature 
of our model although the present errors are too 
large to draw a definite conclusion supporting 
the evidence of / (600) in the B ± ->■ tt+tt~K ± 
decay. Also the preliminary data of the BaBar 
Collaboration [6] below 600 MeV show some en- 
hancement of the tt + tt^K ± events over the back- 
ground. Let us remark that both collaborations 
have not included the a meson in their fits to 
data. Integrating the tttt spectrum between the 



7r + 7r~ threshold and 700 MeV we find the av- 
erage branching ratio for the B + — > aK + and 
B~ — > crK~ decays equal to 3.9 x 10 -6 . A sur- 
plus of events near 0.8 GeV is due to the B ± — > 
p°A ± decay which are not taken into account in 
this model since we concentrate ourselves on the 
S-wavc 77 + tt~ events. 

In figs. 2 and 3 we have shown the averaged 
7r + 7r~ spectra of the two decays: B + — > tt + it~K + 
and B~ — > ir + ir~K~. They have been calcu- 
lated for the charming penguin amplitudes C (m) 
of eq. (6) fitted in ref. [21] by P 1 = (0.068 ± 
0.007) exp[i(1.32 ± 0.10)] and Ff IM = (0.32 ± 
0.14) exp[i(1.0± 0.27)]. Note that the strength of 
the Pi contribution to C(m) is about ten times 
larger than that of Pp /M . If we neglect both of 
them the branching ratios arc smaller by a factor 
of about 4. 

Using these long distance contributions we have 
found a very pronounced direct CP asymmetry in 
the 7r + 7r - spectra. There are many more decays 
of B + into 7r + 7r~A' + than of B~ into -k + ii~ K~ 
(see fig. 3) This asymmetry for the above choice 
of parameters is even higher than the direct CP 
asymmetry recently found in the B° — > K + ir~ 
and B — > K~n + decays [34,35]. The charge 
asymmetry is defined as 



A C p = 



dr(B 




dr(B + ->7r + 7r 


-K+) 




d rriTTTr 


d ra^^ 




dF{B 


--n&ir-K-) 


dT(B+^T&TT 


-K+) 




d m-K-K 


dm^n 





(19) 



If we integrate dT/dm n „ between 0.9 and 1.1 GeV 
then Acp = —0.52 ± 0.12. The errors come from 
the uncertainties of the charming penguin ampli- 
tudes determined in [21]. It would be very use- 
ful to confront this number with a future exper- 
imental determination of this asymmetry in the 
B ± -► /o(980)A ± decays. 

The charge CP violating asymmetry is very 
sensitive to the magnitude and the phase of the 
charming penguin contribution. Using the dif- 
ferent approach of [22] and the fits presented 
in [24] for the B — > Kir decays the value Pi = 
(0.08±0.02) exp[-i(0.6±0.5)] has been obtained, 
while pGiM nas no t been determined. With this 
value of Pi and with pp JM = 0, one obtains 
a good agreement with the HFAG branching frac- 
tion for x = 23.5 GeV -1 . Then the charge asym- 
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mctry is positive and equal to 0.20 ±0.20. We see 
that the different charming penguin amplitudes 
fitted to data arc not in mutual agreement. How- 
ever, the two analyses are based on different data 
sets. In [21] 18 different branching ratios have 
been fitted and in [24] the fit has been performed 
on 8 observables for the B — > Kir decays. 

In the following subsections the predictions for 
the two penguin amplitudes are given without any 
readjustment of the constants \. 

4.2. B° -> 7V+7V-K decays 

In fig. 4 we show the comparison of the model 
predictions for the neutral B decays with the 
BaBar results [4]. Here the experimental back- 
ground contribution is added to the theoreti- 
cal part calculated from eq. (7). The average 
branching ratio for the B° and B decays into 
(ir + ir~)sKg, for the ir + ir~ mass between 0.85 
and 1.1 GcV, equals to 2.93 x 10~ 6 . Twice this 
value compares well with the experimental de- 
termination of the BaBar Collaboration B(B° — ► 
/o(980)AT°) x B(/ o (980) -> ir+ir-) = (6.0 ± 0.9 ± 
0.6±1.2) x 10 -6 . Our curve in fig. 4 is normalized 
to the total number of events multiplied by the 
ratio of 2 x 2.93/6.0 = 0.98. 

The branching ratio diminishes by a factor of 18 
if the charming penguins amplitudes are omitted. 
This was expected due to the near cancellation, 
mentioned in section 2, between the two penguin 
diagram contributions. The absence of tree dia- 
gram in neutral B decays explains the difference 
with the charged B decays where the branching 
ratio drops by a factor of 4 if the charming pen- 
guin terms are not present. 

The direct CP violation asymmetry between 

the decay of B° and B into {ir + -K~)gK s defined 
as 

B (B°^ 7T+ n-K° ) - B(B ^ tt+ n~ K° ) 
_4 = __: : : L (20) 

B (bU 7T+ TT-K° ) +B (B ^ 7T+ 7T- K° ) 

is much smaller than the asymmetry (19) for the 
B ± — > {tt + tt~)sK ± decays. It amounts to 0.01 ± 
0.10 when calculated with the charming penguin 
parameters of [21] in the range between 0.85 
and 1.1 GeV. The reason of its smallness is due 
to the absence of the tree diagram contribution 



40 n — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — i — 1 — 1 — 1 — 1 — r 

35 - 




Figure 4. Effective tt + tt~ mass distribution for 
the B° — > ir + ir~K° decays. The BaBar data and 
the dashed line corresponding to the background, 
are taken from [4]. The solid line is our model 
result. 



(fig. la) for the B° or B decays. 

The BaBar [3] and Belle [8] values for this 
asymmetry, A = 0.24 ± 0.31 ± 0.15 and A = 
—0.39 ± 0.27 ± 0.08, respectively, have large ex- 
perimental errors and agree with our result. The 
HFAG average is -0.14 ± 0.22 [33]. In fig. 5 we 
compare the tt + tt~ spectrum presented by K.-F. 
Chen on behalf of the Belle Collaboration for the 
B° — ► Tr + ir~Kg decays [9] with our calculation 
when the background is subtracted. Our curve is 
normalized to the number of 94 events attributed 
to the /o (980) Kg decay in the range be- 

tween 0.89 and 1.088 GeV [8]. The Belle exper- 
imental determination of the branching ratio for 
the B Q — > tt + tt~ Kg decay is not yet available. As 
in the B ± — * -k + -k^ K^ case, shown in fig. 3, we 
predict some a contribution of the decay ampli- 
tude in the low ir + n~ range below the position 
of p(770) enhancement visible in the data. The 
zero value of our spectrum near 771^^=0.8 GeV, 
also present in fig. 3, comes from = tt/2 (see 
eq. (15)). 

We have also calculated the time-dependent 
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a 




0.6 0.8 
m (GeV) 



Figure 5. Comparison of the Belle data [9] (after 
a background subtraction) with our model (solid 
line) for the B a — > ir + ir~ Kg decays. Dotted ver- 
tical lines delimit a band of the /o(980) events 
used in the curve normalization. 



CP asymmetry in the neutral B° decays into 



A(t) 



dT(B + K% 
dt 



dt 



dT{B ^tt+tt-K°) dr(B°^n + Tr-K°) 



(21) 



dt 



dt 



The time dependence of this asymmetry can be 
approximated as 



A(t) = S sin(Amt) + Acos(Amt), 



(22) 



where Am is the difference between the masses of 
the heavy and light B meson eigenstates. 

In the /o(980) mass range between 0.85 and 
1.1 GeV the asymmetry parameter S equals to 
—0.63 ± 0.09. This result corresponds to the 
charming penguin parameters of [21]. The BaBar 
result is S = -0.95^;^ ± 0.10 [3] and the Belle 
number S = +0.47 ± 0.41 ± 0.08 [8]. The two ex- 
perimental results are not in agreement with each 
other but their experimental errors are large. Our 
result agrees better with the BaBar value. The 
HFAG [33] gives the average S = -0.39 ± 0.26 



which is in agreement with our prediction of 
—0.63. The value of S for the charming penguin 
amplitudes of [24], considered in section 4.1, is 
-0.77. 

4.3. B -> KKK decays 

We have calculated the (K + K~)s spectra and 
the branching ratios for the B + — > (K + K~)$K + 
and B~ — > (K + K~)sK~ decays. For the charm- 
ing penguin amplitudes of [21] one obtains a large 
direct CP violating asymmetry of —0.44 ± 0.12 
in the K + K~ mass range between the threshold 
and 1.1 GeV. The average branching ratio for the 
above mass range equals to (1.8 ± 0.4) x 10~ 6 . 
This value is below the upper limit of 2.9 x 10 -6 
found by Garmash et al. for the branching frac- 
tion B{B+ f a (980)K+, / (980) K+K~) [7]. 
The theoretical (K + K~)s spectrum is flat in 
the range between 1.0 and 1.2 GeV and agrees 
well with the experimental distribution shown in 
fig. 13d of [7]. As in the case of the B^ — > 
tt + tt^K ± decays the asymmetry for the B^ — > 
(K + K~) S K ± process strongly depends on the 
charming penguin amplitude. With the ampli- 
tudes of [24] considered in section 4.1, the asym- 
metry is positive and equal to 0.29 ± 0.21, the 
average branching ratio being (1.7 ± 0.7) x 10~ 6 . 
In these calculations we neglect the symmetrized 
form factors ^2' S ( m KK:) m e 1- (8) assuming that 
for small masses of m K x the fh K -^ masses are 
sufficiently high and that the T2 form factors de- 
crease rapidly with increasing m K -^. 



For the B° 



{K+K-) S K S and 



if 



(K + K~)sKg decays we have obtained very small 
direct CP violating asymmetries A — 0.01±0.10 
and A = 0.001 ± 0.001 for the penguin ampli- 
tudes of [21] and [24], respectively. These num- 
bers agree well with the experimental findings of 
Belle [8] (-0.08 ± 0.12 ± 0.07) and BaBar [36] 
(-0.10 ±0.14 ±0.04). The parameter S, equal to 
—0.64 or —0.77, depending on the set of penguin 
amplitudes, is also in general agreement with the 
results of -0.74 ± 0.27±g || and -0.55 ± 0.22 ± 
0.04 ± 0.11 reported by Belle [8] and BaBar [36], 
respectively. 

Within our model one finds the same asym- 
metries for the B° — * KgKgKg decay as for 
the B° — > (K + K~)sKg process provided that 
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Tabic 1 

Average branching fractions B in units of 10 -6 , asymmetries Acp, A and 5 of our model compared 
to the average values of HFAG [33]. The mass ranges for the B ± -> / (980)X ± and for the 

B° -► / (980)A° decays are (0.9,1.1) GeV and (0.85,1.1) GeV, respectively. The upper limit of the 
[K + K~) s or (KgK^J s effective mass is 1.1 GeV. The model errors come from the uncertainties of the 
charming penguin amplitudes C(m) (eq. (6)) determined in the fits of [21] (model I) or [24] (model II). 
The experimental errors for S in the B° — ► (K + K~)sKg decay are the statistical ones. 





B decay mode 




Average 
HFAG's values 


Model I 
X = 33.5 GeV" 1 


Model II 
X = 23.5 GeV" 1 


B ± 


- /o(980)A±, f -> 7T+TT- 


B 

A C p 


8.49t*;i 
no data 


8.49 (fit) 
-0.52 ±0.12 


8.46 (fit) 
0.20 ±0.20 






B 


6.0 ± 1.6 


5.9 ± 1.6 


5.8 ±2.8 


B° 


-> /o(980)A°, f -> tt+tt- 


A 


-0.14 ±0.22 


0.01 ±0.10 


0.0004 ±0.0010 






S 


-0.39 ±0.26 


-0.63 ±0.09 


-0.77 ±0.0004 




B ± -> (K+K-) s K± 


B 

Acp 


< 2.9 [7] 
no data 


1.8 ±0.4 
-0.44 ±0.12 


1.7 ±0.7 
0.29 ±0.21 






B 


no data 


1.1 ±0.3 


1.2 ±0.5 




B°^(K+K-) S K° S 


A 
S 


-0.09 ±0.10 
-0.55 ±0.22 [36] 
-0.74 ±0.27 [8] 


0.01 ±0.10 
-0.64 ±0.09 


0.001 ±0.001 
-0.77 ±0.0006 






B 


no data 


1.1 ±0.3 


1.2 ±0.5 




B°^(K sK s) s K°s 


A 


0.41 ±0.21 


0.01 ±0.10 


0.001 ±0.001 






S 


-0.26 ±0.34 


-0.64 ±0.09 


-0.77 ±0.0006 



the S-wave is dominant in the production of the 
KgKg pairs and their effective masses arc not 
large. The new data of Belle [37] (A = 0.54 ± 
0.34±0.09, S = 1.26±0.68±0.20) and BaBar [38] 

(A = 0.34±°;28 ± °- 05 > S = -°- 71 -0.32 ± °- 04 ) 

agree for A and disagree for S. Our results for S 
are close to the BaBar value, however the exper- 
imental errors of both collaborations are still too 
large to make a definite conclusion. 

5. Summary and outlook 

We have analyzed the charged and neutral 
three-body B meson decays into the ir + ir~K, 
K + K~K and KgKgKg systems. Our primary 
aim is a construction of decay amplitudes includ- 
ing the Tt + Tr~ final state interactions in a rather 
wide range of effective masses between the tttt 
threshold and 1.2 GeV. The model is based on 
a factorization approximation with some QCD 



corrections and contains the dominant charm- 
ing penguin terms. Using a single amplitude, 
we are able to describe simultaneously the pro- 
duction of two scalar- isoscalar resonances /o(600) 
and Jo (980). The B decay amplitudes to the ttttK 
and KKK states are connected as the coupling 
between the irir and KK channels above 1 GeV 
is incorporated in our model. No adjustable free 
parameters, like arbitrary phase factors between 
contributions of different resonances, are needed. 
We have obtained a good agreement with most of 
the recent BaBar and Belle data for the irir effec- 
tive mass distributions, the branching ratios and 
the time-dependent CP violating asymmetries. 
Our results are summarized in table 1. These 
numbers depend only weakly on the choice of the 
renormalization scale, the corresponding quark 
masses mi, and m s , and the value of F^~^ K (Q). 
The changes are smaller than the errors in the de- 
termination of the charming penguin amplitudes. 
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If we use the charming penguin amplitudes de- 
termined in [21] then the direct CP violation 
asymmetry in the charged B decays to irnK is 
strongly negative (~ —0.5). However, for the 
charming penguin amplitudes taken from [24] this 
asymmetry is positive (~ +0.2). A similar dif- 
ference is found for the B^ — > K + K^K ± reac- 
tions. Future independent measurements of the 
B ± — > tt+tt-K* 1 and B ± — > K + K^K ± asymme- 
tries will be crucial for a decisive test of the phase 
of the long distance contributions. Let us stress 
the importance of charming penguin amplitudes. 
If we omit them then the average branching ratio 
for the B — > ir + Tr~Kg decay is too small by a 
factor of 18 and that of the B ± — > it + it~ mode 
by a factor of 4. Even an ad hoc adjustment of the 
constant x to fit the experimental charged B de- 
cays does not allow one to explain the neutral 
B decays. However, when the charming penguin 
amplitudes are included, we get a good agreement 
for both channels for x values close to the estima- 
tion based on the /o(980) properties. 

The model presented in this paper can be ex- 
tended to larger effective 7T7T mass range, in par- 
ticular to the range where the Jo (1500) is impor- 
tant. One can also include the final state inter- 
actions between one of the pions and the kaon in 
the B or D decays to the irwK system. Especially 
interesting are the n~ K + or n + K~~ subsystems 
where scalar and vector resonances can play an 
important role. 
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